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Abstract. We consider the asymmetric simple exclusion process in one dimension 
with weak asymmetry (WASEP) and 0-1 step initial condition. Our interest are 
the fluctuations of the time-integrated particle current at some prescribed spatial 
location. One expects a crossover from Gaussian to Tracy- Widom distributed fluctu- 
ations. The appropriate crossover scale is an asymmetry of order y 7 ^, times of order 
e~ 2 , and a spatial location of order e _3//2 . For this parameter window we obtain the 
limiting distribution function of the integrated current in terms of an integral over 
the difference of two Fredholm determinants. For large times, on the scale e~ 2 , this 
distribution function converges to the one of Tracy- Widom. 



1 Introduction 



The asymmetric simple exclusion process (ASEP) on the one- dimensional lattice Z 
is a stochastic particle system with at most one particle per site. An ASEP particle 
waits a unit exponentially distributed random time and then jumps to the right with 
probability p, < p < 1, and to the left with probability q = 1 — p. The jump is 
actually carried out only if the destination site is empty. Waiting times and jump 
probabilities are independent. We will set p < |, so predominantly particles move 
towards the left, and always consider 0-1 step initial condition for which the left 
half lattice is empty and the right half lattice is completely filled. We label particles 
from left to right. If x m (t) denotes the position of the m-th particle, m — 1,2,..., 
then x m (0) = m and x m {t) < x m+ i(t) for all t > 0. 

In a celebrated work [I], Johansson investigated for the case p = 0, the totally 
asymmetric simple exclusion process (TASEP), the time- integrated particle current 
across the origin, denoted here by J~(0,t). —J~(0,t) is simply the total number of 
particles which have jumped across the edge (0, 1) up to time t. Johansson proves 
that in the limit t — > oo it holds 

.7(0,t) = -|t + 2- 4 /¥/W (l.l) 

£tw is a Tracy-Widom distributed random variable, which first appeared in the 
context of the large N limit of the largest eigenvalue of a GUE random matrix [2||3] . 



The asymptotics (1.1 ) has recently been extended by Tracy and Widom 4[ 5 to the 



partially asymmetric simple exclusion process (PASEP). In this case J~(0,t) is the 



number of signed jumps across the edge (0, 1) up to time t and, in fact, (1.1) still 
holds with the only modification that t on the left hand side is replaced by t/(q — p). 
In the limit of symmetric jumps, p = q = |, one finds 



J(0, t) S (2tt)- 1/4 (v / 2 - l) Vi t 1/4 eG (1-2) 

for large t, where £g is Gaussian distributed with mean and variance 1. Even finer 
details are established and we refer to the recent study |6] on the large deviations 
for J(0,t). 

With this perspective it is of interest to understand in more detail the crossover 



between the Gaussian central limit theorem (1.2) and the (non-Gaussian) Tracy- 
Widom statistics ( |1.1[ ). The appropriate crossover scale was already identified by 
Bertini and Giacomin [7] . For such a study it is convenient to introduce the dimen- 
sionless scale parameter e, e > and £ <C 1, and to consider the time scale e~ 2 t, 
t = 0(1). The strength of the asymmetry is encoded by the choice 

p=\(l-(3e a ), /3>0,a>0, (1.3) 

which corresponds to a weakly asymmetric simple exclusion process (WASEP). The 
standard WASEP is the particular case a = 1. The macroscopic density profile, 
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p(x,t), is then governed by the dissipative Burgers equation 







ld_ 
2~dx' 



In [8j[9] the associated Gaussian fluctuation theory has been developed which, in 
particular, proves that 



J e (0, e - 2 t) = 



;i.5) 



for small e. The superscript e for J £ should remind that the jump probability 
is adjusted according to (1.3). The variance c(t) 2 can be computed in principle 



from fluctuating hydrodynamics 10 which arrives at an expression in terms of the 



linearization of (1.4) around the time evolved step profile. 

Following [7] the crossover scale is a = |. This crossover scale has also been 
noted in the spectral ga p [ll] of the WASEP generator and in the large deviations 
of the total current 12 13 . We expect that, with the appropriate adjustment of 
constants, (1.5) holds whenever a > |. On the other side for a < | the limiting 
distribution should be Tracy- Widom. In this paper we investigate only the crossover 
regime a = 1/2, for which an added interest comes from the relation to the KPZ 
equation, see (7[[l4] and the discussion in our companion papers 15 , 16 . Hence we 
fix 

p=I(l-/V£), q-p = /3V^, r = ^ = l-2/3 v / i. (1.6) 



On the time scale e 2 t the average time integrated current is of order e 3//2 and a 
typical particle profile increases linearly over the interval [-/3fc~ 3 / 2 , (3te~ 3/2 ] with 
to the left and 1 to the right of this interval. The one-point distribution of the 
time-integrated current will be studied at a general location and not only at the 
origin. 

In Section[5]we summarize our main result, which states that in the rescaled units 
the fluctuations of the integrated current have a size of order t 1//3 with a t-dependent 



amplitude £ t of order 1. The distribution function of £ t is given in (4.29), from which 



it is easily checked that £ t is Tracy- Widom distributed in the limit t — > 00. Thus at 
the crossover scale, one still has the same long time behavior as for the PASEP. 

The analysis heavily relies on the methods developed by Tracy and Widom in El. 
In Section [2] we employ the Ramanujan summation formula, an observation which 
will be instrumental in the asymptotic analysis. The saddle point is discussed in 
Section [3j while in Section [4] we study the /i-integration and convert the contour 
integrations to a Fredholm determinant in L 2 (IR) with a real symmetric kernel. In 
Appendix [C] we argue that one could also perform the \x- integration in the very first 
step, still to arrive at the same result. 
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2 The Tracy and Widom contour integration for- 
mula for the particle's positions 

As discussed in [5], the time-integrated current is directly linked to the motion of 
an ASEP particle with its label properly adjusted. Hence, our focus will be on the 
motion of particles. To ease the comparison we follow closely the notation in |2J, 
which will be referred to merely as TW. We introduce the dimensionless parameter 
a, < a < 1, to label the reference point. At time e~ 2 t the particle index of interest 
equals 

m = a(3te- 3/2 (2.1) 
and the m-th particle is typically at the location 

ci£" 3/2 , a = {-l + 2Va)/3t. (2.2) 

In fact we will have to include a subleading correction to Cie~ 3 ^ 2 of order e~ x I 2 loge. 
The fluctuation scale is 

Cae-i/2 f C2 = a -i/6 (1 _ ^2/3^)1/3 (2 3) 

For later use we also introduce 

c 3 = a-^ii - v^y° /3 m 1/3 , c A = v^/(i - v^) , 

lt = 2/3c 3 C4 = 2/3(/3t) 1 / 3 ( v ^(l - V^)) 2/3 ■ (2-4) 

With these conventions our task will be to study the limit e — > of the WASEP 
distribution function 

Ff(s) = F(x m (e- 2 t) - Cl e- 3 / 2 - (c 2 ht)e- 1 ' 2 \og(2^e) < c 2 se- 1 ' 2 ) . (2.5) 

TW start their analysis with the identity 

lc »k=o' " " ^ 



» oo 

Hx m (t/(q-p))<x)= / n(l-/ir fc )det(l + J(/i))- 

JCa ,._n 



(2.6) 



see TW (25), (27), and Lemma 4. We follow their convention that all contour 
integrals are given a factor l/2wi. Cq is a circle with center at and radius in the 
open interval (r, 1). The operator J(fi) has the kernel J(/x; 77,77') given by 

J(kv,v) = — ri\ ( >w+i — 7 d C 2 - 7 

Jc x Vooiv) w) m+1 C-v 

with fi e C. Here 77,77' are on a circle with center and radius r e (r, 1) and, as a 
linear operator, J(/x) acts on functions on this circle. The integration contour C\ is 
over a circle with center and radius in the interval (1, r/r). (foe is defined through 

V > 00 (v) = (l-Vr x e tM{1 - 1l)) (2.8) 
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and / through 



f(fi,z) = - 



k=— oo 



(2.9) 



By immediate bounds, for \i G C\{0,r n ,n G Z} the function f(fi,z) is analytic 
in the annulus 1 < \z\ < r _1 . We will need its analytic extension, which can be 
deduced from a product formula representation. Following [17] , we set 



i;q)oo = JJ(1 - aq n ) , |g|<l, (a;q) n = (a;q) OQ /(aq n ;q) c 



n=0 



The Ramanujan summation formula, see 17 , Theorem 10.5.1, states 
^ ( a ;q)n x n _ (ax;q) 00 (q/ax;q) 00 (q;q) 00 (b/a;q) 00 



(b; q), 



(x; q) 00 (b/ax; q)^; q)oo{q/a] q) c 



(2.10) 



(2.11) 



provided \q\ < 1 and \b/a\ < \x\ < 1. Setting a — [A, b — fir, x = rz, q — r, one 
easily checks that for 1 < \z\ < r _1 it holds 



(prz; r)oo(l/^; r) oc (r; t^t; r) c 



(rz; r) oc (l/z; t)^; r) oc (r//i; r) 
1 — //z 



oo 

'CO 



(1 _ r n )(1 _ r n) «L (1 - ^ r n)(l - (^)~V>) 

(1 - z)(l - //) Al (1 - ZT «)(1 - z-H-) 11 (1 - /ir")(l - (i-1t») ■ { ■ ' 

Since < r < 1, the right hand side is analytic in z and /i in the domain 
C\{0,r n ,n G Z}. Hence the right hand side of (2.12) is the analytic continua- 
tion of / as defined through (2.9). In both variables / has simple poles at r n , 
neZ. 



3 Saddle point analysis and limit kernel 



We investigate the limit of the kernel ( |2.7 ) of J(/i) as e — > 0. The integrand is 
written as the product of three factors, 



X 



vM){rf) m V'(C-V) 
We study each factor separately and start with Qi. 



x CM) = Qi x Q2 x Q s 



In (2.5) the logarithmic term is switched to the right as 



(3.1) 



'(x m (e- 2 t) - Cl e- 3 / 2 <c 2 (s + \og(2PV~e))s- l/2 ) 



(3.2) 



and thus regarded as a shift of s. First we will ignore this shift, which can be 
included later on because of the uniform error estimates. Hence the parameters in 
the numerator y?oo(C)C m °f Qi are 



m = a/3te~ 3/2 



x = c x e 3/2 + c 2 se 1/2 



(3.3) 



As a consequence, the saddle point analysis is identical to the one of TW with t 
replaced by e~ 3 ^ 2 . The saddle point is given by 



£ = -C4 



see (2.4) for the definition. Setting 

voo(oc m = ^oo(6r^ (c) , 

in a neighborhood of £ = £ one finds, see TW (30), 

^(C) = -(c 3 ) V 3/2 (C - 073 + c 3 *T 1/2 (C - 

+0{ e -*''(c-tf) + 0(e- 1 '\c-Q 2 ). 

The rescaling close to the saddle point corresponds to the substitutions 
V -> C + cJ 1 ^, rj' -tt + ^y/erf , C 6 + c 3 1 V< ■ 

Then 



(3.4) 



(3.5) 



exp[-iC 3 + |(V) 3 + ^(C-V)]- 



(3.6) 



(3.7) 



(3.8) 



For the limit in (3.8), ( G T( and rf G r,,, where consists of the two rays from 
— C3 to — C3 + ooe ±27n//3 , while consists of the two rays from to +ooe ±7r4//3 . 



Close to the saddle point the second factor of (3.1) reads 

c 3 



c 4 ^(C - v) 



(3.9) 



The factor Q 3 needs more work. By ( |3.7[ ) the ratio close to the saddle point 
reads ^ 

4 -l + (c 3 C4)- 1 (V-C)v / i + ^). (3.10) 



£ + c 3 We 



Hence, according to (3.1), /i/(/i, •) has to be evaluated at 1 + ^fez with 

z = (c 3 c 4 ) _1 (V - C) • 



(3.11) 



Correspondingly in (2.12) we substitute z by 1 + yfez and write the product as 
Q 4 x Q 5 x Q 6 . Then 

fif(ji, 1 + V«0 = QzQbQg ■ (3.12) 
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We study the limit of each factor as e — > 0. 
The factor Q 4 reads 

i-mi+^) = i 

To study the limit of the factor Q5 we introduce the g-gamma function, 

T q (x) = j^^(l-q) 1 - x , when|g|<l, (3.14) 



sec 



17 (10.3.3). Setting q x = 1 + yfez one arrives at the identity 

(1 - r n )(l - r n ) 



L (1 - (1 + v^)t™)(1 - (1 + v^^-V") 

= -r T (x)r r (-x)(l + y/iz)-\y/iz)\\ - t)- 2 . (3.15) 

Since q = r = 1 — 2/3y/e, it follows that x = —(z/2/3) + 0(y/e). The g-gamma 
function converges to the gamma function, T, in the limit q — > 1. Hence the limit 
e —7- on the right hand side of (3.15) becomes 

-r(^/2/3)r(-^/2/3)(z/2/3) 2 . (3.16) 

Using 



- zT(-z) = T{-z + 1) , r(*)r(l -z) = , (3.17) 

sin (irz 



yields 



lim fT (l-r")(l-r») _ tts/2/3 

(!-(! + v^^X 1 - (1 + y/ez)- l T n ) sin(7T2/2/3) 



(3.18) 



The factor Q 6 reads 

Q 6 = exp [ £ ( ^ + log 7— , 

n=l r ' 

= e Xp [ E (log (1 - Vi^) + log (1 + ^|z^))] • (3.W) 

71=1 

Since /x G C\R + , the argument of the log lies in C\R_. Hence log is understood as 
the main branch of the logarithm on C\R_. 
Let us define the domain 

V K = \fieC\ sup I m I < 1 , sup |— ^— |<1 + k) (3.20) 

I. 0<j/<1 1 - W) 0<y<l fJ>-y > 



Figure 1: The contour r^uf^. The dots are the poles at (2(3y/e) V and (2/3 y/E) 1 . 



with k > 0. If in (3.19) we assume that fi G P K , then the logarithm can be expanded 
since < r < 1. Hence the exponent [■] in (3.19) reads 



z 



dy- 



dy- 



The error term R(e) is bounded by 
\R(e)\ < c^~e(l/2f3)\z\ 2 



dy 



+ R(e) 



(i - 



2/3 



log(-//) + R(e) 



(3.21) 



+ 



dy 



(v-y)' 



(3.22) 



with some constant c independent of e. According to the definition (3.12), we 
conclude that, uniformly on T> K with an error of order y^loge, 



lim —s/eafia, 1 + \fiz) 



7r/2/3 



sin(7rz/2/3) 



1 (z/2 J 8)log(-M) 



(3.23) 



Since r = 1 — 2j3y/e, the product in (2.6) converges to a non-degenerate limit 
only if fi = 0(y/e). Therefore, as the final step, we substitute \x by 2f3n\fe. Then the 
fi- integration is over a circle with center and radius in the interval (2/3y / e) _1 (r, 1). 
This contour is deformed to U f*. = {//|dist(/z,M+) = 1} n {9fy* < £_1/4 } 
and r^j closes the contour, see Fig. 1, where in brackets we remark that e~ 1//4 and 1 
are taken here only for concreteness. As explained in Appendix A, one can choose 
Eq > such that, for all < e < Eq, fl {/x|3?/x < k} C V k . Also converges 



as e — > to T M = {//|dist(/i, R + 
fi G T M fl {/x|3?/i < ft} it holds 



1}. By uniformity in e one concludes that for 

vr/2/3_ 
/?) 

(3.24) 



1 



/£sm(7rz/2(3) 



x exp [— (log^) + log(2/fyi))] (l + 0(>/i)) 



with z defined in (3.11). 
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To complete the argument, one introduces the logarithmic shift of s and sets 

m = al3te~ 3/2 , x = c^" 3 / 2 + c 2 (s + j^ 1 log(2/3v^))e~ 1/2 . (3.25) 
The rescaled kernel of J(/x) reads then 

J e (/i; T), rf) = J(2/3fx^ £ + c^rjV^, Z + C3 V v^K 1 • (3.26) 



Combining the saddle point asymptotics (3.8) with (3.24), it holds pointwise 



\imJ £ (p;r ] ,r ] ') = I(p; V , V ') (3.27) 

E-S>0 

for fj, G T E , dtp, < K, where the limit kernel is defined by 

V, rf) = J exp [ - f C 3 + HvT + «(C - rf)] ^ 

7T 

sin(7 i " 1 7r(?7 / - C)) 



17 ^(V-OlogC-M)^-!, 



x . - _ 1 - ; — e 7 ' [v-Qiogi-M^-i^ ^ 2g ) 



for // g r M . 

As discussed in Appendix B, the operator J(/i) is trace class for \i G Y il . Hence 
its Fredholm determinant is well-defined. What we would like to show is the validity 
of the limit 

lim det(l + J e U)) = det(l + I (a)) . (3.29) 

e-s>0 

For this we would need the convergence J e (/i) — > as e — > in trace norm. TW 
have to handle the same problem for the limiting case log(— p) = and the sine 
expanded as r y^ 1 n(r]' — £), which formally corresponds to jt — > °°- They use that 
dt(rj' — C) > 0, hence their kernels are bounded. In our case the kernel is singular at 
rf — C — lt n i n G Z. This makes the issue of convergence in trace norm somewhat 
delicate. 



In the following we assume the validity of the limit in (3.29). We also assume 
the exponential bounds 

I det(l + J e {p))\ < ce CoM , I det(l + < ce CoM (3.30) 

for some constant Co < 1. 

4 The //-integration, Fredholm determinant 



Rescaling the //-integration of (2.6) as in (3.26), one obtains 



r 00 1 

F t £ (s)= TT (1 - 2/W^r & ) det(l + J £ (/i))-d/i. (4.1) 

J(2B^)-lCo u-n ^ 
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For fi e C\R_t_ and small e it holds 



- ^r k ) = exp 



fc=0 



d?/-log(l - /jy) 



(4.2) 



Thus only if = (^(a/F) there is a non-degenerate limit and 

oo 

lim Y[ (1 -2f3n^r k ) =e-". 



(4.3) 



fc=0 



Larger values of are exponentially suppressed as exp [— 1 / ^fe ] . We now choose k 
so large that by assumptions (3.30) and by (4.2) the error is small uniformly in e. 
Then 

e- /i det(l + /( / u))-d/i = F 4 ( S ), (4.4) 



limits) 

£^0 



which defines the limiting distribution function i*t(s) still depending on the rescaled 
time parameter t. 

In principle (4.4) is already the final answer, but we still have to transform to a 

us introduce the kernel 



more manageable form. Let us introduce the kernel 

1 



K(ViXuVjXj) 



Q - m 

1 

d - Vi 



it k 



sin( 74 V^- - Q) 



exp [s{Cj-Vj)]F{nj,Cj)- 



(4.5) 



We expand the Fredholm determinant of (4.4). The n-th term of the expansion 
reads 

d/i-e" M / dr] 1 ...dr] n [ d& . . . dC„ 



} det{K(r]i, Q; rjj, Q}i,j=i,...,n ■ (4.6) 



For this expression the /i-integration can be carried out. We set 



w = it 1 ^2ivj - o) • 



(4.7) 



and note that > 0, since for Q e T^, rjj e T v it holds 3?(r/j — Q) > 0. Therefore 
e-^-M/i = 7^7 (e 1 ™ - e- 1 ™)r(w) 

(4.8) 



1 




~~ 27i 


- e 


1 








~~ 2vri 
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In addition 



Q - m 

since ^(Q — r/i) < 0. 

Let us first define the operators with integral kernels 



(4.9) 



df] ^ c d ^sin(7 t ^{q-Q) 



exp [ - |C 3 + W 



+(V - 77Z + 7 f ^ - C) logw ± vy t 1 tt(t] - ()] . 



(4.10) 



We expand the determinant in (4.6) into cycles. The following rearrangement, illus- 
trated for a 3-cycle, is performed for each cycle. We consider the summand e 17TW v w 
of (4.8) using the idendity (4.9) and the definition (4.5) of F. Then the term of 
interest reads 

r 3 

d Vl d V2 d V3 J dCidC 2 dC3n{ e " 7rlfe-0) ^ lfe ^ )F (^'0)} 



x 



POO 

dxi / dx 2 / dx 3 (-l) 3 e Xl(C2 - r?l) e a;2(C3 - r ' 2) e X3te -' ?3) 



dxi / dx 2 / dx 3 (-l) 3 / dr/idCie 



>"7i 1 (»?i-Ci) 1 ,7i ^m-Ci) 



x^(^i,Ci)e 



(Ci^s-^i^i) 



d772dC 2 e i " 7t ~ 1{r?2 " C2 V' 1{m - C2) F(r] 2 , ( 2 )e 



(C2X1-TI2X2) 



x / dr/gdCse^ 7 *" 1 ^ 3 -^^ 7 ' lfe - f3) F(r/ 3 ,C 3 )e 



(C3£2-»?3Z3) 



-l) 3 / dza / dx 2 / dxgtf+fo, x^K+fa, x 2 )K+(x 2 , x 3 ) . (4.11) 



The summand e~ 17rw v w of (4.8) results correspondingly with K£ replaced by K~. 

We thus expand the Fredholm determinant in (4.4), with the n-th term of the 
expansion given in (4.6), and integrate over /1. According to (4.8) this yields the 
two terms corresponding to ±e ±mw v w . Using the identity (4.11) and resumming the 
series results in the difference of two Fredholm determinants. In this difference the 
constant term 1 cancels. Altogether one therefore obtains 

F t (s) = [ e^det(l + /(/i)) -d^ 

= 1 + — / dv-e~ v (det(l - K+) - det(l - K~)) . (4.12) 

Here the determinant is understood in L 2 ([s, 00)) Note that — > as v — > 0, since 
3^(77 — () > 0. Hence for small v the integral in (4.12) is well-defined. 
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Next we want to reexpress the kernels Kf- as Airy- like kernels. Note that, since 
= I\ and T* = T v , one has 

K+(x,yy = K-(x,y). (4.13) 

We also note that, for v — 1, 

K+(x,y) - K-(x,y) = i(27r/ 7t )Ai(x)Ai(y) (4.14) 
by the standard representation of the Airy function, denoted by Ai. Hence 

Kf(x,y) = B t (x,y) ± i( 7 r/ 7t )Ai(x)Ai( 2 /) , (4.15) 
where B t (x,y) is a real kernel and defined through 

B t {x,y)= ( dC / d7 ? (7r/ 7t )cot(7 t - 1 7r(7 ? -C))exp[-|C 3 + ^ 3 + Cy-^] • (4.16) 
Let D = d/dx + d/dy. Then 

[ d( [ d V (r] - C) exp [ - |C 3 + lv 3 + Cv ~ r}x] 

= -D f dC / d?7exp[-|C 3 + ^ 3 + Cz/-^] 

= - J dX5{X)^Ai{x + A)Ai(y + A) . (4.17) 

By functional calculus this identity can be extended to a general class of functions 
of D. To obtain the cotangent, a convenient representation is the power series 

1 °° 1 

TfWfrr 1 ™) = - + 2zJ2 z2 _ (7fw)2 • (4-18) 

The sum converges except for z £ (7 t Z)\{0}. Then 

B t (x,y) = K M (x,y) + J dX'G(0, X')Ai(x + X')Ai(y + A') (4.19) 
with Kai the Airy kernel 

POD 

K Ai (x,y)= dAAi(a; + A)Ai(y + A) (4.20) 
Jo 

and G(A, A') the kernel of the operator 

+ (4.2D 



dX 

n=l 
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By direct computation one verifies that, for a > 0, the operator 2(d/d\)(—(d 2 /d\ 2 ) + 



,2\-l 



has the kernel 



G a (X, A') = -0(A - A')e- a|A - y| 



(4.22) 



where 9(X) = —1 for A < and 9(X) = 1 for A > 0. Summing over n as in (4.21) 
yields 

G(\, A') = -9{X - A , )(e 7 * |A - y| - I)" 1 . (4.23) 

Hence one arrives at 
B t (x,y) (4.24) 

POO 

= K Ai (x, y) + dA(e 7tA - l)" 1 (Ai(x + X)M(y + A) - Ai(x - A) Ai(y - A)) . 



Lemma 4.1. Let B be a trace class operator on the Hilbert space H with scalar 
product (■, ■) and P^ the unnormalized projection along ip GK. Then for a£R 



1 

2i 



( det(l - 5 - iaP^,) - det(l - 5 + iap^,)) 



= -a(det(l - P + - det(l - P)) . 
Proof. Let (1 — B) be invertible. Then 

det(l - B + iaP^) = det(l - P)(l + ia(V», (1 - P)~V)) 



(4.25) 



(4.26) 



Hence the difference on the left side of (4.25) reads 



-adet(l-P)(l + (V, (1-P)~V)-1) = -a(det(l-P + P^)-det(l-£)). (4.27) 
Taking limits on both sides the invertibility condition is removed. □ 



To return to (4.12) we substitute j t u = — logv, du = — (jt v ) 1 df. Then 

K±(x,y) = B t (x + u,y + u)±i(n/ lt )Ai(x + u)Ai(y + u). (4.28) 



Using (4.25) with a = n/jt, one arrives at the final result 



F t (s) = 1 



exp[— e ' 



x 



(det(l - P u (B t - P Ai )P u ) - det(l - P u B t P u ))du . 



(4.29) 



Here P u projects onto [u, oo), Pai has kernel Ai(x)Ai(?/), and the determinants are 
in L 2 (IR). In Appendix C we show that B t is trace class and that the integral in 



(4.29) is well-defined. 
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5 One-point distribution of the particle current 

Let us first state our main result concisely. We consider the WASEP with asymmetry 
/3-s/e and have established that 

\imF(x m (e~ 2 t) - Cl e- 3 / 2 - (c 2 / 7t )^ 1/2 log(2^) < c^ 1 / 2 ) = F t (s) , (5.1) 
where 

m = a/3te- 3/2 , a = (-1 + 2y/a)Pt , c 2 = a- 1/6 (l - ^) 2/3 (f3t) 1/3 , 

7t = 2/3(/3t) 1 / 3 ( v ^(l - v^)) 2/3 , 0<(T<1. (5.2) 



The limit distribution function i*t(s) is defined in (4.29) and its t-dependence is only 
through the parameter jt. 

We reexpress our findings in terms of time-integrated currents and set 

J7" e (j,t) = ft of signed jumps across the bond + 1) up to time t, 

where the superscript e reminds on the ^-dependence of the asymmetry. The trans- 
formation from particle position to integrated current is discussed in jH], which can 
simply be followed (their quantity X equals — J). We also allow for a shift of the 
reference point by e~ 1 x, x = 0(1). Then, for \y\ < 1, we define the (x, independent 
family of distribution functions 

Ff x>t) (s) = F(2^J%[yf3te- 3 / 2 + xe- 1 \,e~H) + |(1 - \y\) 2 fi 2 te~ l 

-(1 - \y\)We~ 1 ' 2 + (x 2 /2t) - log(2/3v^) < lt s) (5.3) 

with |_-J denoting integer part. In the new parameters 

lt = 2- 1 / 3 (l-y 2 ) 2 / 3 /3 4 / 3 t 1 / 3 . (5.4) 

foFF xt Js) is the properly scaled probability distribution for the statistics of the 
current at location [y/3te~ 3 ^ 2 + xe~ l \ integrated over the time span [0,e~ 2 i]. The 
limit (5.1) can then be restated as 

hmF^ t) (s) = F t (s) . (5.5) 



Let £t be a random variable with distribution function F t . Then (5.3) translates 
to the one-point statistics of the time-integrated current as 

2(3^eJ e ([y/3te~ 3/2 + xe' l \ , e~ 2 t) 

= -|(1 - W^te- 1 + (1 - {yDWe- 1 / 2 - (x 2 /2i) + \og(2^) + 7*6 , (5.6) 

valid in the limit of small e. 
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It is instructive to write down the probability density for £ t , which is given by 



j roc 

Pt(s) = ^F t (s) = J 7 *e 7t(s - u) exp [ - e*<-% t (u)du , (5.7) 

where 

g t (u) = det(l - P u (B t - P Ai )P u ) ~ det(l - P u B t P u ) . (5.8) 



1 — exp[— e 7 ' s ] is the Gumbel distribution function and the first factor in (5.7) is 
the Gumbel probability density. Since lim s _ 5 ._ 00 Ft(s) = 0, it follows that the second 
factor of the convolution is normalized as 



g t (u)du = l. (5.9) 
However, numerical computations clearly indicate that the definite sign of pt is 



regained only after smearing with the Gumbel density, see 15 



As discussed in the Introduction, in the long time limit the integrated current 
i s ex pected to be Tracy- Widom distributed, which cor responds to taking t — > oo in 
(5.6). Since jt ~ t 1 ^ 3 , the Gumbel distribution in (5.7) tends to 8(s — u) and B t of 
(4.24) tends to Kx\ as t — > oo. Hence ^ has the probability density 

Poo(s) = goo(s) = F TW (s)u(s) (5.10) 

with 

F TW (s) = det(l - P s K Ai P s ) , u(s) = (P,Ai, (1 - P a K M P a )- 1 P a M) . (5.11) 
It follows from the identities in [3] that 

(log F TW )' = u. (5.12) 

Hence = F^ w and 

lim F t (s) = F TW (s) , (5.13) 

t— >oo 

confirming the conventional expectation. 



6 Conclusions 

We have studied the crossover asymptotics of the WASEP in an e-dependent pa- 
rameter window. Somewhat more physically, our result can be rephrased through 
appropriately adjusted time scales. One considers the asymmetry q — p which is 
assumed to be small but fixed, (q — p)^ 1 defines the time unit. For short times the 
statistics of the integrated current is approximately Gaussian. In an intermediate 
time window one should observe the statistics corresponding to F t , which for long 
times approaches the Tracy- Widom distribution -Frw As q — p is increased the 
crossover window should shrink and might become not discernible at all. 
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In this paper the focus is on the derivation of F t (s) from the WASEP scaling 
limit. Properties of this family of distribution functions and their relation to the 



KPZ equation will be discussed in the companion papers 15, 16 



Remark: After submitting this article, G. Amir, I. Corwin, and J. Quastel posted 



their paper 18 in the arXiv, in which independently they establish the limit (5.1) 



and the formula (5.7) for the probability density pt- As here, their starting point is 



the PASEP Tracy- Widom contour integration formula. 

Acknowledgements. We are grateful to Michael Prahofer for many illuminating 
discussions and to Sylvain Prolhac for helping with the combinatorics in Appendix 
C. H. S. thanks Jeremy Quastel for emphasizing the importance of the crossover 
WASEP. This work is supported by a DFG grant. In addition T.S. acknowledges 
the support from KAKENHI (9740044) and H.S. from Math-for-Industry of Kyushu 
University. 



A Appendix: The set V K 

For k > we define V K = n V {2) , with 



W = U sup \-J^—\<l\,V^ = U\ sup |-JU<i + l6 l (A.l) 

0<y<l ' 1 - PV '0<j/<1>-2/' 

and T K = {/x|dist(/x, R+) = l,9fyx < k}. We claim that one can choose eo > such 
that for all < e < Eq it holds 

2/V^r K cD K . (A.2) 

For n G 2(3y/eT K one has |9fyt| < 2(3(1 + k)^/e, < 2(3 y/e. Clearly, for 

sufficiently small e our assertion holds for the set T>^ 1 '. 
To discuss we define 

i{y)=y{{a-y) 2 + b 2 Y l, \ 0<y<l. (A.3) 

(i) For a < 0, the maximum of £ is at y$ — 1 and it holds that £ < 1. (ii) For a > 
and a 2 + b 2 < a the maximum of £ is at y = (a 2 + b 2 )/a and I < |6| _1 (a 2 + 6 2 ) 1 / 2 . For 
jtt G 2(3 ^fe T K we set /i = 2(3y/e(a +ib). If a < 0, we are in case (i) and our assertion 
holds for this part of the contour. If a > 0, then a = 2(3^a, a < k, b = 2(3\feb, 
\b\ = 1. For e sufficiently small it holds 2(3\fe{a 2 + b 2 ) < a and we are in case (ii). 
Hence £ is bounded by b^ 1 (a 2 + fe 2 ) 1 / 2 < (k 2 + l) 1 / 2 < 1 + k, which implies that also 

(2) 

this part of the contour is contained in T> K . 



B Appendix: Bounds on the kernels I(fi),K v 



± 



The (^-integration in J(/i) from (3.28) has a simple pole whenever 



rj' — ( = 7 t m , mGZ. (B.l) 
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Thus the (^-integration is singular for a discrete set of points rf G T v . The integration 
along looks locally like 



du(u — z) 



with sufficiently small a > 0. For z 7^ and 
uniformly bounded and extends by continuity to z 



(B.2) 



by c e 



-%l%-%f 



I < a/2, the integral (B.2) is 
0. Thus \I(n',T},rf)\ is bounded 



for sufficiently small 5 > 0. Hence I(fx) is trace class. 



A similar discussion applies to the kernels K^(x,y) from (4.10). The singularities 
of csc(7T 7 f ~ 1 (?7 — £)) are simple isolated poles only, hence integrable, the integrations 



in (4.10) are well-defined, and the kernel is bounded by c^e 
small 5 > 0. 



-%I 3 p-%T 



for sufficiently 



C Appendix: //-integration first 



We describe an alternative route to arrive at (4.12). The starting point is the expan- 
sion of the Fredholm determinant (2.6), which is a 2n-fold integral over (1, . . . ,(n 
and 771 , ... , r) n . Each integrand has factors corresponding to Qi, Q 2 , Q3 as in (3.10). 



For each of them we use the saddle point approximation as discussed in Section [3j 
There then remains an n-fold product of functions / with an argument scaled as in 



(3.10), (3.11). Finally the //-integration has to be carried out, see (2.6) 



The new idea here is to first integrate over /x and then take the limit e — > with 
the arguments \fe close to the saddle. Thereby one circumvents the discussion at 
the end of Section [3] With a more careful variant, possibly, one could control the 
error bounds. 

To lighten the notation, we use instead of r the conventional symbol q = 1 — \fe. 



For 1 < 



, n, we define 



f 1 - 

H n (z u ...,z n )= / dfji-g q (n) TTW(/x, Zj)} 

J Co ^ j= i 



„ n 

-IT E / d^g^H^-q-^iz^}, (C.l) 



with 



9q(f>) = Hi 1 - /V 



(C.2) 



k=0 



We split the sum over 
Thereby 



into n-tuples with no double points and the rest. 



H n = H}*' + HI 



(C.3) 



Only Hn is discussed. We checked that the remainder term H™ st tends to for 
e — > for n = 2, 3, 4, 5. Unfortunately the combinatorial structure becomes involved 
and we did not try to work out the extension to general n. 
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Hn^ has only the simple poles q e , £ G Z. C encloses those poles with £ < 0. 
Therefore the /^-integration can be carried out with the result 

n 

h£\ Zi , ...,z n ) = (-1)" w $>({4 < o})?"^^^) 

(€i,...A)ez™ »=i 

xf[<HQ-* ~ Q- ei )fl(*i) ei ■ (C4) 

Here £^ means no double points and means that the i-ih factor is omitted 
from the product. We now substitute £j by £j + £i for j ^ i. Then 



oo n 

H^z u ...,z n ) = (-ir{Y i g 9 ^)ll(z j ) t ) 



£>0 j=l 
n n 

i=1 (< 1 ,..,4..,<„)e(z\o)"- 1 i= l 
= HW(z 1 ,...,z n )Hl*^(z 1 ,...,z n ). (C.5) 

We first study the issue of analytic extension. Since g q (q e ) = exp[— cq e ], Hn^ ,T 
is analytic on all of C ra . For H^'^ it suffices to consider the case i = n. From 
the restriction of no double points, we conclude that we have a sum over products, 
where each factor is of the form 

oo 

^(1 - q- m )- k z m = {(-!)*(! - q m Y\q k z) m + (1 - g m )-V m } (C.6) 

m^O m=l 

oo oo 

^i,...A=o m=i 

oo 

^i,...,4=o 

for = 1,2, The last expression is the analytic extension to C \ {0} except for 

the poles at z = q e , £ G Z. 
We summarize the result in 

Proposition C.l. Let V a = {z G C | 1 < \z\ < q" 1 } and let V q = {z G Z | z ^ 

0,z 7^ q e ,£ G Z}. Then Hn \ as defined on T>® n , extends to an analytic function on 

<? ' 

Close to the saddle z 3 - — 1 + \fewj. For 1 + \fswj G T> q we define 

fj(*) _ rr(*),r rr(*),S //-i ~\ 

JJ n,e -"ji,e JJ n,e ) v^- 1 ) 

Hif r ( Wl , ...,w n ) = {-l) n {^e) l - w H^ T {l + y/iw u . . • , 1 + V£w n ) , (C.8) 
H^(w u ...,w n ) = (V^y-'H^il + y^i, • • • , 1 + >fav n ) , (C9) 
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with the shorthand w = Y^j=i w r We a ^ so introduce the limit functions defined on 
V® n , with V = C \ Z, 



H^(wi, . . . , w n ) = 7T _1 sin(7ra>) ( J^J 7T _1 Sin(7TWj)J 



(CIO) 
(C.11) 



Theorem C.2. Pointwise on T>q it holds 



lim fl"W ,r = #5 , lim #!*}' s = fZ? 



e-s>0 



e— >0 



(C.12) 



The proof is divided into several parts. We start with g q . 
Lemma C.3. It holds 



lim(v^)- +1 Yl 9 9 tf) ( + ^0) ~ l = r < 



w 



(C.13) 



3=1 



Proof. The product in flC.13| equals 1 + y/ew + 0(e 2 ). Then, by [17], Eq. 10.3.1, 

oo oo 



£W 



£=0 

r i (g; g)oo 
(y;g)oo 



, y = (1 + . (C.14) 



The g-Gamma function is defined by 

T q (x) = (l-q) 



i-x (?;s)c 



(? x ;?)c 



|g| < 1 . 



(C.15) 



Setting x = w, one thus arrives at 

oo 

(V^)- w+1 Y,g 9 tf)(i + V^w)- e = (i + v^r'r 



(C.16) 



^=1 



which implies the limit (C.13). 



□ 



Lemma 3 establishes the left part of (C.12) 

Lemma C.4. The following limits hold, 
(i) For k = 1 



iTo^E { - (1 - (1 + V~ew 1 )q e+i y\ 1 + V~ew 1 )q e+1 + ((1 + sTewi) ~ <?TV} 



£=0 



7T COt TTWi . 



(C.17) 
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(ii) For k>2 



lim( v ^) fc V {(-l) fc (l - (1 + y/Ew 1 )q k+tl+ - +t "y (1 + V~e Wl )q k+ ^ + - + ^ 



£ lv ..,4=0 



+((l + ^i)- g 4+ - +4 )"g^ + - +4 } 
1 



(1 + (-l) fc ) 



(fc-l)! 



dim^fe" 



(*-l)! 



B(k). (C.18) 



For even k the coefficients -B(fc) are related to the Bernoulli numbers Bk by 
B(k) = k- 1 (27r) k (-l) 1+ ^B k . 

Proof, ad (i): Separating the I = term of the second summand and expanding 
inside the curly bracket one obtains 



limVe[- (1 - (1 + y/ewi)) - ^ {(* + V^i)(l - y/ef) (l - (1 + y/ewx) 

x(l - ^Tel)y l - (1 - ^i)(l - ^)(1 - (1 - v^iXl - V~sl)Y 1 }] 

= 71 COt TlWi . (C.19) 

ad (ii): For k > 2 the sum approximates the Riemann integral 



/•oo /*oo k k 

/ dux... / dn fc r[e-^(l-r[e-^)" 1 (l + (-l) fc ). 
Note that the limit does not depend on Wi, . . . , Wk, in contrast to = 1. 



(C.20) 



□ 



With Lemma 3 and 4 we have identified the limit of i?n*e , in principle. It 
remains for each n to rearrange the sum such that the expression (C.ll) results. 

To handle the constraint of the sum in (C.5), let us denote by J n _i the set 
{1, . . . , n — 1} of labeled vertices. An unoriented edge with endpoints i ^ j, is 
denoted by b and 5& stands for S^. Let £ n -\ be the set of all edges of I n _i. Then 
the constraint of the sum in ( |C.5[ ) can be written as 



n c 1 -*)- 



(C.21) 



By expanding the product one has to sum over the set Q n -i of all undirected graphs 
over I n -\. For g E Q n -\ the weight, w(g), in this sum results from (— l)W ed s es ) and 
from the limits (C.17) and (C.18). A given graph g decomposes I n -i into disjoint 
clusters. In view of (C.ll) and (C.17) we keep the number of clusters of size 1 
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fixed (they yield a product of cotangents) and sum over all other clusters, by (C.18) 
necessarily of even size. Thereby we arrive at the following counting problem. 

Given is the set I m of to vertices, m even, and set Q m of undirected graphs over 
J m . For a given graph g e Q m , I m decomposes into r clusters Ci,. . . ,C r of size 
| Cj- 1 = rrij, Y^j=i m j = m - The weight w(g) is defined by 

r 

W (g) = (_!)»(«*«) JJS( mi )(( mi _l)!)-l. (C.22) 

J'=l 

Note that B{m) = for odd m. Since we do not want to allow for rrij = 1, only for 
the next lemma we set 

5(1) =0. (C.23) 

Lemma C.5. For even m it holds 

J2 W ^) = — r(-l) m/2 - (C24) 

g&Qm 

Proof. The weight w(g) induces a weight of the clusters as 

r 

w({C 1 ,...,C r }) = l[w(C j ). (C.25) 

To compute w(Cj) we note that 

J2 w{g) = (k - l)!(-l) fc - 1 ((^ - l)!)- 1 ^) • (C.26) 

g£5fc, single cluster 

The prefactor can be verified firstly because it holds for k = 2. Now assume it is 
valid for general k > 2. Then, adding an extra vertex, k + 1, it can be connected in 
distinct ways to the cluster of size k and the number of edges is thereby always 
increased by 1. Hence 

w(C J ) = (-ir- 1 B(m J ). (C.27) 



We introduce a generating function, /(A), for the left hand side of (C.24) by 

fw= E srE»w 

even m=2 g(zGm 
oo oo 



E E^ E *(E-*-o 

mi,...,m r >2 
r 

(-irj]^( 



to! r. 

even m=2 r=l m\,...,m r >2 j=l 



m 

mi . . . m r 



[m j/ 



3=1 



oo oo 1 

EwH) r (E/W- (C28) 



r! ^— ' A;! 

r=l fc=2 
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Using (C.18) the sum over k reads 

POO 

/ du(u(e u - l))- 1 2(cosh Xu - 1) = log(vr A/ sin vr A) . (C.29) 
Jo 



x | 



Hence 

/(A) = (7rA)- 1 sin(7rA) -1. (C.30) 

Taylor expanding / confirms the claim. □ 

Proof of Theorem C.2. We return to with the constraint in the summation 
written as in (C.21). For g G Gn-i we decompose I n -\ = I*_ x U I^-i- In-i consists 
of one point clusters {i}, i G I^-i, an d I^-i consists of clusters of size > 2, compare 
with the notation above (C.22). By (C.17) the one point cluster {i} carries the 
weight 

hi = 71 COt(7TU>j) . (C.31) 

Either set could be empty and \In-i\ is even. With this convention the weight of g 
is given by 



w(g) 



r 

-ip d ^( n /^n^KXK-i)!)- 1 }, (c.32) 



where the second product refers to the clusters of see (C.22). If = 0, then 
the first product equals 1 and correspondingly for I^-i- Then 



H^(w u ...,w n ) = ^/ n _i(wi, . . . ...,w n ) 



i=l 



with 



(C.33) 



(C.34) 



f n - 1 (wi,...,w n -i)= ^ w (9)- 

It is convenient to introduce the symmetrizer S n . If g is some function on C 
1 < m < n, then S n g is defined by 



(S n g)(w 1 ,...,w n ) = — ^2g(w n{1) ,... 



i ^7r(m) ) 7 



(C.35) 



where the sum is over all permutations n of (1, . . . ,n). Let us also set \I^-x\ — j, 
j — 0, 2, . . . , [n], where [n] = n — 2 for even n and [n] — n — 1 for odd n. In the sum 
(C.34) we fix the set and sum over all graphs for the set I^_ v By Lemma C.5 
this yields 

(C.36) 



( n '^jTT-'f- 1 



U'/2 
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Performing the sum over all subsets I^-i, one obtains 



/ - 1 \ 

f n „ 1 (w ll ...,w k _ 1 ) = £ (-l^V^ + l)" 1 n-l-j JSn-i^i---^-!-,)- 

j=0,j even ' 

(C.37) 



Carrying out the sum (|C.33|), it then follows that 
H^(w u ...,w r , 



») = E (-l) i7 V ( n H _. \ S n {h X ■ ■ ■ h n +j) 

•i — l"l rein N / 



J=0, even 



TlWq 



= 7r 1 sin (7r ifj) y Y\ 71 1 sin( 
The last term is a mere rewriting of the identity 



mx>)= e (-i) (fc - i)/2 e n^n^- 



odd fc>l 
k<n 



Ac{l,...,n} ieA i£A c 
\A\=k 



(C.38) establishes the second assertion of Theorem C.2. 



(C.38) 



(C.39) 



□ 



If the remainder term in (C.3) is ignored, then by Theorem C.2 the scaled H n 
converges to 



n n 

H^H^ = r^vr" 1 sin (tt E w i) ( II 71-1 sin( 



7TWj 



(C.40) 



j=i i=i 



This agrees with (4.6) upon performing the //-integration by using the first identity 



of (4.8) and collecting the factors from the saddle point. 



D Appendix: Trace class property 

Proposition D.l: The operator P s B t P s is trace class. The functions appearing in 



(4-29) and (5.1) are absolutely integrable in u. 

Remark. We have no direct proof that pt(s) > and pt(s)ds = 1. 



Proof. All operators will be defined on L 2 ([s, oo)). 
We have 

B t = K M + C t , (D.l) 

where for simplicity we set 74 = 1. In general, if B = A1A2, then (tr|5|) 2 < 
(tr^Ai)(tr^A 2 ), see pi, Section VI. 6. For the Airy kernel we write 



/oo 
dAAi(x + A - s)M(y + A - s) . (D.2) 
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Hence = A\ with Ai(x, y) = Ai(x + y — s) and 

rco poo 

tr|X Ai | < dx dy\Ai(x + y) | 2 . (D.3) 

J s J s 

For the operator C t we write 

poo 

C t (x,y)= / dA(e A — (Ai(x + \)Ai(y + A) — Ai(x — \)Ai(y — A)) 

/oo 
dA(e A - s - ly^Aiix + A - s) - Ai(x - A + s))Ai(y + A - s) 

/oo 
dAe -(A-s)/2 Ai ^ _ A + s ) ( Ai ( y + A - s) - Ai(y - A + s)) 

x(e A ~ s -l)- 1 e (A - s)/2 . (D.4) 

Hence C t = A 2 A X + A 3 A 4 . 
By definition 

poo />oo 

tr|A 2 | 2 = / dx dy(e y -iy 2 (Ai(x + y)-Al(x-y)) 2 . (D.5) 

We split the y-integration into the intervals [0, c], [c, oo). In [0, c] we Taylor expand 
Ai(x + y) — Ai(x — y) in y and choose c so small that the first order dominates. We 
then determine C\ such that (e y — 1)~ 2 < c 1 e^ y on [c, oo). Using Schwarz inequality 
yields the estimate 

/oo . pc 

dx( y dyyV - l)- 2 Ai'(x) 2 

dyAi(x + yf + 2 J dye~ y A\{x - yfj (D.6) 
for a suitable choice of the constant c 2 . By definition 

/oo /»oo 
dx y dye- y Ai{x - yf (D.7) 

and, repeating the argument for tr | ^4 2 1 2 , 

/oo poo 
dx y d^e" - l)-V(Ai(x + y)- Ai(x - y)) 2 

dxi^J dyy 2 (e y -l)- l e y Ai'(x) 2 

/oo poo 
dyAi(x + yf + 2 J dye~ y Ai(x - y fj . (D.8) 

Using that |Ai(x)| ~ exp[— x 3//2 ], |Ai'(a;)| ~ exp[— x 3//2 ] for x — > oo and |Ai(x)| ~ 
|:r|~ 1//4 , |Ai'(a;)| ~ Ix] 1 ^ for x — > — oo, all integrals are bounded with the asymptotics 
exp[— s 3,/2 ] for s — > oo and |s| 3 / 2 for s — > — oo. 
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Let us consider the integral (5.7) for pt(s). The first factor decays as e~ u for 
u — y oo and as exp[— e'"'] for u — > — oo. For the second factor we use the inequality 
det(l + B)\ < exp[tr|_B|], see [20], Section XIIL17. From our previous estimates 



on the trace norm of B u gt(u) is bounded by c for u — > oo and as cexp[|w| 3//2 ] as 



u — > — oo, which establishes integrability. For the integral (4.29) defining F t (s) one 



uses that for large u the determinants behave as 1 + 0(exp[— u 3 ^ 2 ]). □ 
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